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1.  INTRODUCTION 


\ 


Advanced  design  of  aerodynamic  and  hydrodynamic  vehicles  often 
depends  critically  upon  an  understanding  of  transition  sensi- 
tivity to  many  phenomena.  Nosetlps  on  ballistic  reentry 
vehicles,  for  exam{.ie,  must  be  designed  to  withstand  transi- 
tion destabilizing  effects  of  surface  roughness,  ablation, 
surface  cooling  (in  the  presence  of  roughness),  and  entropy 
gradient.  As  a second  example,  current  advanced  hydrodynamic 
vehicle  design  procedures  make  advantageous  use  of  the  strong 
stabilizing  effects  of  favorable  pressure  gradient  and  sur- 
face heating  to  generate  laminar-flow  vehicles.  Linear  sta- 
bility theory  is  the  most  popular  theoretical  tool  for  guiding 
design  of  many  vehicles  for  which  transition  location  influ- 
ences the  design  configuration.  While  linear  stability  theory 
has  proven  irrelevant  for  reentry  vehicle  nosetlps  (because  of 
the  occurrence  of  roughness-induced  transition  bypass),  a 
great  deal  of  success  with  stability  theory  has  been  enjoyed 
for  hydrodynamic  vehicles,  most  notably  by  Wazzan  and  Smith. ^ 

For  transition  triggered  by  small  amplitude  disturbances, 
linear  stability  theory  provides  a nearly  exact  solution  to 
the  Navler-Stokes  equations.  There  is  li^’tle  doubt  that  sta- 
bility theory's  Tollmien-Schllchting  waves  exist  and  play  an 
important  role  in  the  initial  stages  of  transition.  Because 
the  end  state  of  the  transition  process  is  a (highly  nonlinear) 
turbulent  flow,  linear  stability  theory  breaks  down  at  some 
point  between  that  of  the  Initiation  of  Tollmien-Schlichting 
waves  and  the  transition  point  (defined,  for  example,  as  the 
point  where  skin  friction  achieves  a minimum) . In  other 
words,  linear  stability  theory  is  inapplicable  in  the  post- 
critical  stages  of  transition  and  therefore  has  no  natural  way 

^ of  specifying  the  actual  transition  point. 
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Undaunted  by  this  limitation  on  stability  theory's  applicabil- 
ity.  Smith  and  van  Ingen  simultaneously  (with  no  knowledge 
of  the  other's  activities)  devised  the  well-known  e®  method. 

As  is  so  often  the  case  with  clever  approximations,  the  e® 
method  has  yielded  accurate  predictions  for  flows  well  beyond 
the  original  data  base.  However,  the  record  of  success  has 
been  blemished  somewhat  by  inaccurate  predictions  including, 
surprisingly,  applications  such  as  transition  sensitivity  to 
freestream  turbulence  with  varying  spectral  content. 

Presumably,  either  a nonlinear  stability  theory  or  an  exact 
Navier-Stokes  solution  method  is  needed  to  rigorously  bridge 
the  gap  between  the  initiation  of  Tollraien-Schlichting  waves 
and  the  transition  point.  Because  neither  of  these  two 
approaches  has  been  developed  to  the  point  of  being  practic- 
able for  engineering  design,  the  designer  must  depend  upon 
existing  correlations  and/or  approximate  methods  such  as  the 
e®  method. 

Recently,  a new  approximate  transition-prediction  method  has 
been  devised  which  shows  great  promise  for  engineering  design, 
viz,  the  turbulence-model  transition-prediction  method. 

The  method  is  based  on  the  conventional  long-time  averaged 
equations  of  motion.  Nonlinear  correlation  terms  such  as  the 
Reynolds  stress  are  approximated  in  a manner  similar  to  that 
used  in  standard  closure  schemes  for  turbulent  flow  modeling. 

As  will  be  explained  in  greater  detail  in  Section  2,  turbu- 
lence-model equations  on  the  one  hand  are  expected  to  apply 
in  the  latter  stages  of  transition.  On  the  other  hand,  the 
time-averaging  process  removes  explicit  appearance  of  Tollmlen- 
Schlichtlng  waves.  Consequently,  the  turbulence-model  transi- 
tion-prediction method  is  reasonly  well  founded  only  near  the 
end  of  transition. 
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As  with  the  llnear-stabllity/e’  method,  turbulence-model 

equations  have  provided  accurate  transition  predictions  well 

beyond  the  original  data  base.  In  fact,  the  equations  formu- 

7-9 

lated  by  Wilcox  and  Chambers  have  yielded  accurate  results 
for  virtually  all  of  the  applications  made  to  date.  However, 
some  of  this  success  has  been  achieved  with  an  adjustable 
parameter,  viz,  the  freestream  turbulence  level.  Additionally, 
the  method  has  no  natural  way  of  representing  spectral  effects. 

In  summary,  because  correlations  generally  are  limited  to  a 
restricted  data  base,  1 inear-stabillty/e®  and  turbulence-model 
methods  are  the  only  two  comprehensive  transition  theories 
which  are  of  practical  utility  for  the  aerodynamic/hydrodynamic 
vehicle  designer.  The  former  is  theoretically  sound  only  dur- 
ing the  Initial  stages  of  transition  while  the  latter  is  well 
founded  only  near  the  end  of  transition.  Hence  we  reasonably 
may  speculate  that  a synthesis  of  these  two  methods  will  yield 
a transition-prediction  theory  which  is  fundamentally  more 
sound  than  either  theory  standing  alone. 

The  primary  objective  of  this  project  has  been  to  synthesize 
1 inear'-stabl  lity  theory  and  the  turbulence-model  transition- 
pr*- diet  ion  method.  In  so  doing,  the  most  immediate  result 
of  the  proposed  research,  if  successful,  would  be  development 
of  a physically  sound  alternative  to  the  empirical  Smith-van 
Inge::  e’  procedure.  Aside  from  the  advantage  of  having  a more 
fundamental  method  for  simulating  the  post-critical  stages 
of  transition,  the  proposed  synthesis  would  obviate  the  costly 
eigen-solutions  needed  to  compute  amplification  ratios  beyond 
the  critical  Reynolds  number.  Rather,  a straightforward 
boundary-layer  marching  computation  would  be  all  that  is  need- 
ed in  solving  the  turbulence/transition  model  equations  from 
the  critical  Reynolds  number  up  to  the  transition  Reynolds 
number . 
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As  will  be  shown  in  the  following  sections,  an  important 
first  step  has  been  taken  toward  accomplishing  the  proposed 
synthesis.  Results  for  the  Blasius  boundary  layer  show  that 
the  most  important  (in  the  context  of  transitional  boundary 
layers)  closure  coefficient  appearing  in  the  turbulence-model 
equations  approaches  a universal  limiting  value  for  amplifica- 
tion ratios  in  excess  of  about  e**  . 
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2.  ANALYSIS 

This  section  first  reviews  the  turbulence-model  equations 
which  form  the  basis  of  the  study.  Then,  the  key  closure 
coefficient  controlling  transition-prediction  accuracy,  A, 
is  cast  in  terms  of  linear-stability  variables.  Finally, 
the  coefficient  X is  computed  for  the  Blasius  boundary  layer 
over  the  entire  stability  spectrum  and  is  found  to  approach 
a universal  limiting  form  for  amplification  ratios  in  excess 
of  e**  . 

2.1  REVIEW  OP  THE  TURBULENCE-MODEL  EQUATIONS 


2.1.1  The  Model  Equations 

For  arbitrary  incompressible  flows,  the  model  equations 
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Turbulent  Dissipation  Hate 


pu 


3u)‘ 


J 3x 


■ !“  ■^2Sjj  - [b  + 2o(||-)‘>t  P»' 


3x 


J L 


(u  + ope) 


3a) 

3x 


il 


(5) 


where  Xj  is  position  vector.  Time-averaged  (mean)  velocity 
Is  denoted  by  Uj  while  h,  p,  p and  y are  mean  enthalpy, 
pressure,  density  and  molecular  viscosity;  and 

are  the  shear  stress  tensor  and  heat-flux  vector.  The 
quantity  Sj^j  = 1/2  (3Uj^/3xj  + 3Uj/3x^)  Is  the  mean  strain- 
rate  tensor.  The  turbulent  mixing  energy,  e,  and  the  turbu- 
lent dissipation  rate,  o),  are  needed  to  define  the  eddy 
dlffuslvlty,  e,  which  Is  given  by  the  following  equation: 


e = e/u) 


(6) 


In  order  to  close  this  system  of  equations,  the  stress  ten- 
sor and  heat  flux  vectors  must  be  specified.  In  the  Wilcox- 
Tr'ac  1 tur’bulence  model  the  stress  tensor  Is  assumed  to  have 
Its  tn’lnclpal  axes  alllgned  with  those  of  the  mean  sti’aln- 
rate  tt'nsor  so  that  we  write' 


- 


ij 


(7) 


w 


Altjo,  appealing  to  the  classical  analogy  between  heat  and 
momentum  transfer,  we  write  the  heat  flux  vector  as 


- / v*  + p£  N 

"VPrj^  Pr^j 


(8) 


whre  Prj^  and  Pr^  are  laminar  and  turbulent  Prandtl  numbers. 
The  quantity  H Is  the  turbulent  length  scale  defined  as 


Z 


e^/u) 


(9) 


The  turbulent  Prandtl  number,  Pr^,  and  the  closure  coefficients 
a,a*,S,B*,<7,a*  appearing  in  Equations  (4)  and  (5)  are 


R = J- 
^ 20 


fi*  = — 2— 
100 


a = 


O.  . -- 


Pr  = — 
T 9 


(10) 


a = -y[l  - (1-A)  exp  (-Re^/2)] 
cx»  = ^[1  - (1-X)  exp  (-2Re^)] 


where  Re^p  is  the  turbulent  Reynolds  number  defined  by 


Re^  = pe^£/y 


(11) 


As  will  be  discussed  in  the  next  subsection,  transition  pre- 
dictions are  most  sensitive  to  the  closure  coefficient  X 
appearing  in  the  last  two  of  Equations  (10).  Detailed  study 
of  the  viscous  sublayer  of  a turbulent  boundary  layer^’^^  indi- 
cates that  accurate  sublayer  properties  can  be  simulated 
with  (see  Figure  1) 


X 


1 

11 


(12) 
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Dissipation  -<-H — ► Production 
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Utilrif’;  Kquatiot^s  (1)  through  (12)  (In  some  cases,  similar 
vi^rsions  of  these  equations),  a wide  range  of  compressible 
and  incompressible  turbulent  flows  iiave  been  computed, 
Including  complicating  effects  such  as  boundary-layer  separa- 
tion,^^ streamline  curvature, sui’face  roughness^^  and 

mass  injection.  With  no  additional  modifications,  tiie  model 
accurately  simulates  many  aspects  of  incompressible  flat-plate 
boundary  layer  (FPBL)  transition.  Hence,  at  least  for 

the  relatively  simple  FPBL,  this  model  for  fully  developed 
turbulent  flows  provides  an  accurate  description  of  mean  flow 
properties  during  transition.  And  this,  despite  the  fact  that 
the  long-time  averaging  process  leading  to  Equations  (1) 
through  (5)  masks  the  presence  of  any  wave-like  instability! 
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Before  proceeding  to  further  discussion  of  why  the  model  does 
as  well  as  It  does  and  of  the  value  of  X for  more  complicated 
t'lows,  it  is  instructive  to  present  the  physical  defitAitions 
of  the  turbulence  parameters  e and  o)  and  appropriate  surface 
boundary  conditions  for  each. 

7 1 H 

As  argued  by  Wilcox  and  Chambers,  ' ^ tlie  tui’bulent  mixing 
energy  is  proportional  to  the  kinetic  energy  attending,  the 
fluctuation  of  fluid  particles  normal  to  the  plane  of  shear. 
Letting  v'  denote  the  fluctuating  velocity  component  normal 
to  tiie  shear  plane  (under  the  boundary-layer  approximations, 
shear  planes  are  parallel  to  the  x-z  plane),  the  turbulent 
mixing  energy  is  given  by 


e 


(11) 


'Hie  pipYslcal  meaning  of  u)  lias  also  l)ecn  discussed  by  WMcox  and 
7 15 

Chambers.  ’ For  I ncompress lb le  boundary  layers,  comparison 
of  the  limiting  forms  of  tlie  model  equations  and  the  exact 
Heynolds  stress  equation  very  close  to  a solid  boundary 
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q 


showti  that  0)  Is  tht;  rate  at  which  e Is  dissipated  into  heat, 
mean  kinetic  energy  and  other  fluctuation  modes;  a suit- 
able definition  of  u)  is 

„ , iv  <OvV3yr> 

6*  <v''>  (I**) 

Using  the  definitions  given  in  Equations  (13)  and  (1^), 
boundary  conditions  appropriate  for  smooth  and  rough  surfaces 

7 

have  been  devised  by  Wilcox  and  Chambers  and  by  Wilcox  and 
Traci.  Letting  k denote  (peak-to-valley ) roughness  height, 
detailed  study  of  the  viscous  sublayer  shows  that  (in  the 
limit  of  small  roughnesses): 

(15) 

(16) 


2.1.2  Fundamental  Considerations 

The  only  feature  of  transitional  boundary  layers  which  has 
been  explicitly  introduced  into  the  model  is  the  following. 
Concurrent  with  our  analysis  of  the  viscous  sublayer,  the 
value  of  \ has  been  established  by  demanding  that  the  linear- 
stability  minimum-critical  Reynolds  number,  Re^  , for  the 
Blasius  boundary  layer  match  the  corresponding  model-equation 
neutral-stability  Reynolds  number,  Re~.  The  latter  Is  defined 
as  the  Reynolds  number  based  on  the  plate  length  beyond  which 
turbulent  energy  production,  a*l3u/3ylpe,  exceeds  turbulent 
energy  dissipation,  0*p(i)e.  Using  the  Blasius  profile  and 
noting  that  m-ZOv/Qy^  [see  Equation  (16)],  the  neutral 
stability  Reynolds  number  is 
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e = 0 


at  y = 0 


^ y , smooth  surfaces 
p y* 

V 

4320  — ^ , rough  surfaces 


as 


Demanding  that  = Re~  yields  the  value  of  A given  in 

i'Aluation  (12). 

The  model  equations  reasonably  can  be  expected  to  apply  to 
transitional  flows  which  are  insensitive  to  spectral  effects. 
That  is,  the  various  constants  in  the  model  equations  are 
essentially  correlation  coefficients  which  have  been  inte- 
grated over  the  turbulent  spectrum.  Hence,  if  the  stability 
diagram  shows  that  a wide  range  of  wave  numbers,  a,  undergo 
amplification,  the  spectrum  will  more  closely  resemble  a 
fully-turbulent  spectrum  than  if  only  a small  range  of  wave 
numbers  are  unstable.  For  example,  the  stability  diagram  for 
a boundary  layer  subjected  to  a pressure  gradient  is  shown  in 
Figure  2.  For  adverse  pi’essure  gradient,  a finite  range  of 
wave  numbers  are  unstable  at  all  Reynolds  numbers  in  excess 
of  Rev  (note  that  6*  is  displacement  thickness).  On  the 
basis  of  the  discussion  above,  the  model  would  be  expected 
to  accur’ately  predict  the  destabilizing  effect  of  adverse 
pressure  gradient.  in  contrast,  the  stability  diagram  becomes 
thinner  with  Increasing  favorable  pressure  gradient  so  that 
spectral  effects  become  Increasingly  important,  particularly 
for  small  freestream  turbulence  Intensity,  T’ , which  yields 
transition  at  large  values  of  the  model  hence  would  be 

expected  to  fare  poorly  for  transitional  boundary  layers 
with  favorable  gradients  (and  small  freestream  disturbances). 
The  model  behaves  Just  as  the  above  discussion  indicates. 
Kxcellent  agreement  between  theory  and  experiment  is  obtained 
for  adverse  gradients  while,  for  low  freestream  turbulence 
intensities,  the  model  fails  to  predict  tiie  strong  stabilizing 
effect  of  favorable  gradients. 


To  remove  this  deficiency,  Wllcox^^  has  extended  the  re- 
quirement 

Re-  = Re  (l8) 

* 

to  include  favorable  pressure  gradients,  suction  and  surface 
heating  effects.  For  small  freestream  turbulence  intensity. 
Equation  (12)  is  replaced  by 

(\5/2 

r(A)  (19) 

wher'e  A is  given  by 


and  the  quantity  0 is  momentum  thickness.  The  function  f(A) 
has  been  determined  by  Invoking  Equation  (l8).  Figure  3 
presents  results  based  on  linear-stability  predictions  for 
the  Pohlhausen  profiles  and  for  the  asymptotic  suction  pro- 
file;^*^  a good  fit  to  the  data  indicates  the  variation  of 
the  function  f(A)  with  A is  hence 

f(A)  = ^ + ^exp  r-iJOA  H(A)]  (21) 

where  H(A)  is  the  Heaviside  stepf unction . 

Tn  essence,  Equations  (19)-(21)  represent  a correlation  of 
llneai’-stablllty-predlcted  minimum-critical  Reynolds  numbers. 
We  have  thus  implicitly  built  in  some  of  the  wave-instability 
phenomena  which  were  lost  through  the  long-time  averaging 
procedure.  Two  key  points  must  be  amplified  regarding  this 
last  point.  E'irst , the  coefficient  A is  not  Invariant  for 
transitional  flows;  rather,  it  is  sensitive  to  the  specti-um 
of  unstable  frequencies.  Thus,  it  is  through  the  precise 
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value  of  X that  Tollmlen-Schllohtlng  waves  implicitly  appeal’ 

111  the  model.  Because  the  only  transi tion-speci Tic  modifica- 
tions to  the  modt'l  equat.lons  ai’e  for  the  coefficient  X 
Ll''.quat  1 ons  (19)-(21)1,  It  Is  lienee  unsur-prl  s I ng  that  ultimately 
we  have  chosen  to  rely  upon  linear-stability  theory  to  set  its 
value.  oecond , and  of  central  import-ance  to  the  present  j 

program , the  coefficient  X also  controls  the  rate  at  which  dis- 
turbances are  amplified  beyond  the  critical  point.  The  model's 
ability  to  accurately  predict  actual  ti’ansitlon  point  location 
indicates  that  once  the  critical  point  is  coi’rectly  established, 
th('  model  provides  an  accurate  simulation  of  the  post -crit  1 cal 

* stages  of  transition  with  no  further  appeal  to  stability  theory. 

diven  this  Insiglit,  the  whole  concept  of  using  turbulence-model 
eciuatlons  to  describe  transition  can  be  cast  in  a different 
light.  On  the  one  hand,  time-averaging  conceals  many  impor- 
tant physical  aspects  of  transition  mechanisms,  particularly 
during  the  early  linear-amplification  phase.  We  ai’e  thus  obli- 
gated to  put  some  of  the  physics  back  into  the  equations  which 
was  lost  through  the  t ime-avei’aging  process;  ergo,  the  modifi- 
cation to  X.  On  tlie  othei’  hand,  assuming  the  lat.ter  phases  of 
t.  ranr.  1 1 1 ot;  1, 1 . e . , close  to  the  transition  point)  to  be  very 
rapid,  the  time-averaging  pi’oeess  becomes  a more  plausible 
concept  as  the  flow  more  nearly  resembles  a tui’bulent  flow, 
t interest  Itigly , the  modification  to  X has  little  effect  on  the 
latter  st.ages  of  transition.)  The  turbulence-model  transition- 
prediction  approach  thus  has  its  strongest  foundation  in  the 
latter  stages  of  transition,  pi’eclsely  the  regime  where  coti- 

j vnetlonal  linear-stability  methods  are  not  well  founded. 

I Thus,  turbulence-model  equations,  might  most  properly  be  viewed 

as  a plausible  alteiuiatlve  to  the  e®  metiiod.  That  is,  a llneai’ 

I stability  conquitation  could  be  performed  up  to,  and  perhaps 

a bit  beyond,  the  minimum-critical  Reynolds  number.  Results 

I 


I 


of  the  stability  computation  would  define  Then,  rather 

than  continuing  to  solve  eigenvalue  problems  to  determine 
amplification  factors  up  to  the  e*  amplification  point,  the 
model  equations  could  be  used  to  predict  transition  location. 

As  noted  above,  what  we  are  currently  doing  is  using  a cor- 
relation of  linear-stability  minimum-critical  Reynolds  numbers 
to  fix  the  value  of  A. 

Note  that  there  is  a fundamental  premise  underlying  the  notion 
that  a computation  with  the  model  equations  can  replace  the 
conventional  e®  procedure.  The  premise  is  that  A must  be 
constant  (or  at  least  neai'ly  constant)  for  varying  Reynolds 
number  and  also  be  either  constant  or  assume  a universal 
variation  through  the  boundary  layer.  If  A is  either  constant 
or  assumes  a universal  variation  through  the  layer,  there  is  a 
rational  procedure  for  setting  its  value  in  terms  of  the  com- 
puted eigenfunctions.  Likewise,  if  A is  nearly  constant  with 
Reynolds  number,  there  may  be  a rational  procedure  for  choos- 
ing an  appropriate  "average"  value  such  as,  for  example,  the 
value  at  the  point  where  the  linear-stability  solution  has 
been  amplified  to  e"  times  its  Initial  value.  In  other  words, 
if  linear-stability  theory  indicates  that  A either  (a)  does 
not  vary  substantially  with  Reynolds  number  and  with  distance 
through  the  boundary  layer  or  (b)  varies  in  a universal  manner, 
the  closure  approximation  involving  A will  be  proven  for  very 
low  Reynolds  number. 

2.2  THEORETICAL  FORMULATION 

As  noted  above,  the  overall  objective  of  this  research  has 
been  to  synthesise  linear-stability  theory  with  the  turbulence- 
model  transition-prediction  method.  In  the  proposed  synthesis, 
stability  theory  is  used  to  locate  the  critical  Reynolds  number, 

H . Then,  based  on  the  stability  solution  at  R (or  perhaps 

C '' 

at  a slightly  larger  Reynolds  number),  initial  conditions  are 
established  for  a turbulence-model  computation  which  proceeds 
ftn^m  the  poitU  at  which  the  stability  solution  is  valid,  up  to 

i 

I 
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tht'  tfariisltlon  point.  To  obtain  a well-posed  Initial  value 
pfot)l<‘ni  foi-  thi'  l.urbii  1 1'nce-model  computation,  Initial  con- 
ditions are  nec.-ded  foi*  three  quantities,  viz,  e,  w and 
In  tills  subsection,  these  three  quantities  are  specified 
In  terms  of  linear-stability  eigenfunctions. 


Confining  analysis  to  incompressible  flows.  Equations  (13) 
and  (l4)  provide  an  obvious  way  of  defining  e and  oj  profiles 
in  terms  of  the  stability  theory  eigenfunctions.  To  deter- 
mine the  initial  value  of  X we  must  examine  the  equation  for 
e.  For  incompressible  boundary  layers,  the  equation  for  e 
simplifies  to 


where  u and  w are  the  mean  velocity  components  in  the  stream- 
wise,  X,  and  lateral,  s,  directions.  Note  that  in  writing 
Equation  (22)  conventional  boundary-layer  approximations  have 
been  made  and,  in  addition,  the  flow  has  been  assumed  paral- 
lel so  that  V = 0.  The  closure  coefficient  o*  assumes  a value 
of  1/2  while,  in  the  limit  of  small  Reynolds  number  based  on 
turbulence  intensity  and  scale  [see  Equation  (10)J  , a*  simpli- 


fies to 


= 10  ^ 


We  now  proceed  to  derive  a relation  between  the  closure  coef- 
ficient X and  the  stability  solution.  Noting  the  definition  of 
«,  we  can  derive  an  exact  equation  for  its  evolution  by  taking 
the  v'  moment  of  the  v-momentum  equation.  The  following  equa- 
tion, subject  to  the  same  approximations  used  in  writing 
Equation  (22),  results  for  e. 


•3x 


2 p 3y 


V < 


/3v' 

" 3e 

^ ■ . 3 J 

V9y 

) 9y 

L'ly  ■ 

(23) 


Comparison  of  Equations  (22)  and  (23)  shows  that  we  have  made 
the  following  two  closure  approximations: 
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(24) 


(25) 


With  some  manipulation,  the  closure  coefficients  X and  o* 
can  be  expressed  as  follows  . 


j* 


v<v*3>  <Ov'/9y)^> 

<v'  9<v'  ''>/3y 


(26) 


X 


10 

3 


2<fi  - v< 
p 9y 


^ 

<v’  > 


9u/9y I 


2 

> 


(27) 
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Equation  (2?)  Is  the  desired  relation  which  can  be  used  along 
with  Equations  (13)  and  (l4)  to  define  Initial  conditions  for  a 
turbulence-model  computation.  Potentially,  Equation  (26)  pro- 
vides a further  check  on  the  validity  of  the  closure  approxima- 
tions. However,  as  will  be  shown  below,  Equation  (26)  pre- 
dicts 0*  = 0 as  a consequence  of  the  assumed  linearity  of  the 
solut Ion . 

In  the  linear-stability  solution,  the  velocity  and  pressure 
fluctuations  v'  and  p'  are  written  as 

v'(x,y,z,t)  = U^4)(y)exp[l(ax  + Be  - wt)]  (28) 

p'(x,y,z,t)  = pU^TT  (y  )exp  [l  (ax  + (?z  - (it)]  (29) 
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where  t denotes  time,  is  freestream  velocity,  S and  § 
are  wave  numbers,  (j3  is  frequency,  and  the  functions  (})(y) 
and  TT(y)  are  the  complex  amplitude  functions  of  the  disturb- 
ance flow  variables  v'  and  p'. 

To  evaluate  the  time-averaged  quantities  appearing  in  Equations 
(13),  (1^),  (26)  and  (27)  we  use  the  following  definition: 


H}j>  = 


OoT-voo 


t+T 

^ / i|'(x,y,z, 

m 


t ' )dt ' 


An  immediate  consequence  of  the  assumed  linearity  of  the 
solution  is  that  correlations  of  odd  order  vanish,  e.g.. 


:v ' > = <v ' ^>  = 


■’5>  = ...  = 


Thus,  as  noted  above.  Equation  (26)  implies  that  o*  = 0. 
Although  o*  is  postulated  to  be  1/2  by  Wilcox  and  Traci,  this 
discrepancy  is  of  little  consequence  as  the  diffusion  terms 
in  the  e equation  play  an  insignificant  role  in  transitional 
flows . 

Working  with  the  real  parts  of  the  linear-stability  solution 
functions,  performing  all  time-averages  indicated  in  Equations 
(13),  (1^)  and  (27),  and  denoting  local  Reynolds  number  by  R 
yields  the  following: 


9 . .2. 


(d(t>^/dy)^  + (d<}i^/dy)^ 


+ 


X = - 


1_0 


r>((|)^,ilTi^/(i,V  + (|)j(iTTj/dy)  - H ^ [(il(|)^,/<i,v)‘  + (d(|)j/d.v)‘ ] 

+ <t.‘^)|9U/3y| 


wh<M'('  y 1:'.  now  d 1 mon:’.  I on  1 ('Of.  dlr.tanco  fi'otr.  the  f.url’ncc?,  U = , 

and  oubf.cf  1 ptf.  r and  1 d('not('  foal  an  Imafjnafy  part.  Kquat;  1 otio 
(^;>)  _ (1^1)  ai'o  t ho  dof.lfod  folat  lonr.  whloli  can  bo  uf.od  to  dol’lno 
Init  ial  condlblotif.  for  a t.urbu  1 ('neo-modol  comput  at.  1 on  In  l,»'ftno 
ol'  a p,lvon  llni'ar-f.lab  1 1 Ity  oolut.lon. 


■!  KVAl, HATTON  OK  TIIK  OI.OllUHK  OOK, KKICTKNT  X 

llavlnp  di'flv('d  Kquat  lonf,  abovi'?,  Lhe  next  ot.cp  If.  t.o  per- 

Toi’m  a 1 1 noar-f.t  ab  1 1 1 t.y  oompnt.at  1 on  and  t o ('xamlne  t.ho  feoultant 
pfot’llof.,  pant  loulafly  t.ho  X pi-ol’lli'.  The  mof.t.  aonslblo  ataftliu*, 
P('lnt  la  to  l'0[:':ln  with  the  Hlaf.lna  boundary  layer.  A larpo  number 
or  f.tal'lllty  eotnput  at  1 onf.  have  boon  perfoi’iiu'd  with  the  Mack^"^ 
fdablllty  prot-tram;  all  comput, at.  I ons  have  beeti  done  with  the  spatial 
amp  1 1 1’ 1 oat  Ion  t hooi'y  iiptlon.  both  Hc'ynolds  numbei’  and  frecpioncy 
Irivo  l'(M>n  v.'irled  In  order  to  determine  t.he  variation  of'  X t hrouy;h- 
out  t h<'  Hoyno  Idf.-numl'er/rri'quenoy  plane. 

lO'lor  t('  ('xamlnlnf,  ri'f.ults  of  tlu'  computations,  It.  is  Ins-truot  Ivi 
ti)  roc.-ill  t.hat.  the  basic  prc'inlst'  of  I’ormlnt-’;  a syntheslf-  ot'  linear 
f.t.ablllty  and  turbu  1 ('iice-model  methocis  Is  that  two  condlt.ions  I'o 
sat  1 s (‘ l('d . Tlu'  first  condit  ion  Is  th.at  the  closure  cot^ri'lclent 
\ b('  I'lther  relat  ively  weak  f'unct  loti  oT  dls.tance  1 hrouph  t.he 
t'oundary  layer  or,  at  It'asd.  , assume'  .•)  similar  variation  f’ot'  all 
unstable  I’ri'quenc  1 (?s ; this,  ccuidltlon  valldat.es  t.tu?  clof.uri'  ap- 
pi'e-ix  1 mat  1 e'n  ih'i’lned  In  Kquat  ton  (,'”■>).  The  su'cond  condlt  Ion  1 . ■ 
that.  \ v.ai'.v  .-.lowly  with  Iv'.ynolels  numbf'i*;  this,  e'e'ndltlon  If.  r.oid- 
t'll  (at  t ('  v.i  1 1 d;vt ('  K.qu.-it  Ion  (Oh)  and  (bl  to  ('llmlnate'  .-.oltil  Ion 
f-i'Uf.lt  Ivlty  to  Inlt  l;il  I'e-'iid  1 1 1 onf. . Note'  also  t h.-it  we-  heipi'  te' 
fillet  the'  ;i\e'raKe'  value'  elf  X t.ei  be-  re''.-!  f.onali  1 .V  e'lof.e;  to  the  pe'f.tulal- 
e'  l hlaf.  luf.  V.-llUe',  V 1 f.  , X = 1/11. 


,’0 
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Figure  chows  computed  A profiles  at  nine  Reynolds  numbei'c  cor’- 

respondlng  to  one  stable  case,  one  neutrally  stable  case,  and 

seven  unstable  cases  corresponding  to  amplification  from  t!ie 

n 

neutral  case  by  factors  of  e with  values  of  n ranging  from  0 
to  10;  for  all  nine  cases  the  frequency  is  g,ivt?n  by 

Fr  = mv/U^  = O-IO”^  (S'")) 

W(^  ai’e  thus  following  the  evolution  of  A for  a constant-frequency 
disturbance  initiated  at  a Reynolds  number  upstream  of  the  neutral 
point  corresponding  to  the  frequency  given  iri  Equation  (35 1 ■ 

As  shown,  although  A varies  rapidly  with  y near  the  surface,  all 
nine  curves  display  approximately  the  same  variation.  However, 
above  a vlaue  of  n = y /U^/vx  of  about  ?,  the  various  A profiles 
vary  rapidly  with  n and  do  so  in  dissimilar  manner  for  amplifica- 
tion  ratios  up  to  e . For  example,  at  the  lowest  Reynolds 
number  (for  which  the  solution  is  stable),  A is  negative  above 
n - 2.  As  we  move  to  ttie  neutral  point  we  find  that  A vanishes 
foT'  values  of  n In  excess  of  2.5.  Then  as  Reynolds  number  in- 
creases, A varies  more  and  more  rapidly  with  n and  asymptotes 

to  a single  curve  for  amplification  ratios  in  excess  of  between 

S , 4 

O'  and  e . 

Figure  5 shows  similar  cur‘ves  for  a dimensionless  frequency  Fr 
given  by 

Fr  = 3- i0"5  (36) 

Again,  the  curves  collapse  to  a single  curve  for  amplification 

l| 

I’atios  in  excess  of  e . 

Computations  have  been  performed  for  frequencies  covering  '.he 
entire  stability  diagram.  For  each  fi’equency  considered,  tiie 
computed  A profiles  always  asymptote  to  a universal  pr-oflle 
beyond  the  e point  with  a subtle  qualification.  That  ir,  r-c- 
ferrltig  to  Flgui'e  6,  as  Reynolds  number  increases,  the  uppei' 
brancii  of  the  stability  diagram  eventually  is 


reaciied  and  we 


Figure  6.  Stability  dlagi’am  for  the  Blaslus  boundary  layer 
depicting  a constant  frequency  trajectory;  n de- 
notes solution  amplification  of 


again  enter  a stable  region.  As  we  approach  this  upper  neutral 
point,  the  X profiles  begin  to  fall  back  to  those  typical  of 
low  Reynolds  numbers. 

The  i-a[)ld  var-iation  of  X near  0 = 0 results  from  a breakdown  in 
the  basic  closure  approximat  ions  near  the  sui’face.  That  is, 
tlie  produc.t  ion  term  in  the  <v’^>  equation,  <-p—  to 

zero  quadratlcally  with  distance  from  the  surface  so  tliat , in 
terms  of  n , 

<-~  as  n^O  (37) 

p oy 

By  contrast,  the  modelc'd  production  term  for  He,p->0  behaves  as 

X ~ e ^ X n"  as  n-*-0  (38) 

Consequently,  close  to  the  surface  we  ultimately  have 

X ~ as  tr*-0  (39) 

This  modeling  .shortcoming  i ;>  of  little  consequence  as  dissipation 
exceeds,  production  near  n=0.  Consequently,  for  the  remaindei’  of 
this  discu.asion  our  focus  will  be  upon  the  region  between  n = l 
and  the  outer  edge  of  the  boundai*y  layer,  n=5. 

Klgure  7 shows  computed  X profiles  for  several  frequencies  and 
amp  1 li’icat i on  ratios.  As  shown,  all  of  the  computed  X profiles 
cluster  about  the  approximate  profile  defined  by 


X = .009  3 + .0015exp  [-|^(n-l'>] 


Hence,  one  of  the  basic  I’equi  rements  for  estaM  ishing  a ‘'.ahuiiuice 
modei/1  Inear-stability  synthesis  appears  to  be  satisl'ied,  v'z, 
the  X pr  '’lie  appears  to  appr’oacli  a universal  (1.e.,  independent, 
of  frequency)  limiting  pi'oflle  as  amplification  ratio  inoroasos. 
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Figure  7.  Profiles  of  the  closure  coefficient  X for  various 
frequencies  and  amplification  ratios. 


To  dotor-mlne  the  rate  of  approach  to  the  asymptotic  prol’lle. 

It  Is  convenient  to  examine  the  average  value  of  X defined  as 


fol  low:: 


X dn 


KIgur’e  8 shows  X as  a function  of  displacement  thickness 
Reynolds  number,  note  that  Re^*  denotes  the  neutral- 

stability  value  of  Re^*  for  a given  frequency.  As  shown,  fei’ 

the  higher  frequencies  X approaches  its  asymptotic  value  most 

_S 

rapidly.  For  Fi-  = 0.5‘10  , the  lowest  frequency  at  wliich 

computations  have  been  done,  the  approach  to  the  asymptotic 
value  is  the  least  rapid.  A key  feature  of  all  the  curves  is 
that  their  asymptotic  values  lie  between  X = .066  and  X = .O83, 
as  compai*ed  to  the  postulated  turbulence-model  value  of  .OQl. 

Furtlier  examination  of  the  X variation  with  Reynolds  number  show: 

that  for  the  higher-  frequencies,  the  peak  value  Is  achieved  at 

14 

an  ampl  1 ficat ion  ratio  of  about  e while,  for-  the  lowest  fi’e- 

quency,  X is  about  half  Its  asymptotic  value  at  this  i*atio. 

Figure  9 presents  the  vai-iation  of  X witli  frequency  for’  amplifl- 
4 

cation  ratio  e , including  a correlation  of  the  computed  values; 
the  correlation  Is: 


|l  - exp  [-|^(10*^Fr)]| 


The  asymptotic  values  of  X ai’e  also  shown  for  reference.  Two 

key  conclusiorrs  oarr  be  drawn  from  the  observed  variation  of  X. 

4 - 

First,  beyond  an  amplification  ratio  of  e , X varies  slowly  with 

Reynolds  number’.  Second,  X is  only  weakly  fi’oquency  dependen* 

_5 

foi’  frequencies  in  excess  of  Fr*  = 1.0  • 10  . These  two  points  . 

lotrd  further  cr’edence  to  the  second  of  tdie  two  basic  hypotheses 
underlying  a synthesis  of  turbulence-model  and  1 Inear-stab  11  It. y 
theories . 
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^ DincUll'iVION 


Results  presented  In  Section  2.3  show  that  for  the  Blasius 
boundary  layer  the  closure  coefficient  X approaches  a universal 

4 

limiting  form  for  amplification  ratio  in  excess  of  about  e 

while,  for’  smaller  ampl  1 f icat ion  ratios,  A appears  to  be  sti’ong- 

ly  fr’equency  dependent.  This  observation  is  consistent  with  the 

notlon  that  a uniformly  valid  theory  can  be  developed  by  using 

4 

linear  stability  theory  up  to  the  e amplification  point  and  by 

then  using  the  turbulence-model  equations  from  this  point  up  to 

transition.  On  the  one  hand,  using  stability  theory  up  to  am- 

4 

pilficatlon  ratios  of  e is  quite  reasonable  as  nonlinear’  terms 

almost  certainly  are  negligible  in  this  regime.  On  the  other 

hand,  the  weak  dependence  of  A on  frequency  for  amplification 

4 

ratio  greater  than  e and  the  universal  limiting  form  are  con- 
sistent with  the  overall  notion  of  using  long-term  averaging. 

Thus,  for  the  Blasius  boundary  layer,  we  have  a suitable  de- 
finition for  the  post-critical  stages  of  transition,  viz,  the 

stages  beyond  the  point  at  which  a boundary-layer  disturbance 

4 

has  been  amplified  to  e times  its  Initial  value. 

In  summary,  although  results  presented  represent  only  a first 
cut  at  accomplishing  a synthesis  of  turbulence-model  and  linear- 
stability  theories,  encouraging  progress  has  been  made.  Further 
research  in  two  specific  areas  is  needed  to  complete  the  synthesis. 
The  first  area  needing  further  investigation  is  evaluation  of  A 
for  boundary  layers  with  pressure  gradient,  suction  and  surface 
heat  transfer.  In  so  doing,  we  can  determine  whether  a profile 
such  as  that  defined  in  Equation  (40)  applies  to  all  boundary 
layers  or  if  straightforward  generalizations  can  be  made  to  de- 
vise a universally  applicable  A profile.  The  second  area  is  use 
of  Equation  (40)  in  a tui’bulence-model  computation  to  determine 

4 

solution  sensitivity  to  (a)  the  point  of  initiation  (e.g.,  e 
amplification  point)  and  (b)  freestream  turbulence. 
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HoiUllts  of  this  study  suggest  that  when  this  research  has  been 

done  our  ultimate  goal  of  establishing  a physically  sound 

q 

alternative  to  the  empirical  Smith-van  Ingen  e'  procedure  can 

be  realized.  j 
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